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Abstract 

We obtain a non-linear generalization of the relativistic diffusion. We 
discuss diffusion equations whose non-linearity is a consequence of quan- 
tum statisticsWe show that the assumptions of the relativistic invariance 
and an interpretation of the solution as a probability distribution sub- 
stantially restrict the class of admissible non-linear diffusion equations. 
We consider relativistic invariant as well as covariant frame-dependent 
diffusion equations with a drift. In the latter case we show that there can 
exist stationary solutions of the diffusion equation besides the equilibrium 
solution corresponding to the quantum or Tsallis distributions. We define 
the relative entropy as a function of the diffusion probability and prove 
that it is monotonically decreasing in time when the diffusion tends to 
the equilibrium. We discuss its relation to the thermodynamic behavior 
of diffusing particles . 

1 Introduction 

An interaction of a particle with a many particle system is too complex to be 
described by a microscopic Hamiltonian theory. Some approximations for the 
interaction are unavoidable. If the medium is passive then an approximation of 
a linear master equation is sufficient [1] . In the quantum case there is a restric- 
tion on such an approximation (even if the interaction is weak) resulting from 
quantum statistics. The probability of a transition to a quantum state depends 
on the occupation number of this state. Such a dependence implies a non-linear 
master equation [2] . The collision terms in the Boltzmann equation [3] as well as 
interactions in an equation for the Wigner function [4] in quantum field theory 
also lead to non-linear integro-differential equations. Under the approximation 
of two-body interactions the non-linearities are quadratic. The approximation 
of binary collisions together with a non-relativistic approximation to the proba- 
bility distribution of the electron gas is applied to electron-photon interactions 
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[5] [6] leading to the Kompaneets diffusion equation for the photon distribution. 
The simplicity of the Kompaneets equation inspired us for a search of relativistic 
partial differential equations of the parabolic type. 

In this paper we discuss diffusion equations from the point of view of the gen- 
eral principles of relativistic invariance. First, we consider relativistic invariant 
non- linear diffusion equations in the phase space (x,p). In general, differential 
equations invariant under a group of transformations L have solutions which 
transform in a non-trivial way under L. We consider rotation invariant station- 
ary solutions. They transform covariantly under the Lorentz group. After a 
transformation the solution depends on the velocity of the reference frame 
(the boost). We may treat the frame velocity as an additional physical variable 
and discuss relativistic equations for densities which depend on the phase space 
variables (x,p) and on w. In such a case we obtain a wider class of non-linear 
diffusion equations similar to the Kompaneets equation. 

The linear relativistic diffusions have been studied for a long time (see the 
reviews [7] [8]). It has been shown a long time ago that Markovian relativistic dif- 
fusion in the configuration space does not exist [9] [10]. The relativistic diffusion 
(preserving the mass-shell) of a massive particle on the phase space has been 
defined and studied by Schay [11] and Dudley [12]. A general class of diffusions 
in the phase space has been discussed in [13]. An analog of Dudley's diffusion in 
general relativity has been investigated in [14]. We have studied diffusions with 
a friction of massive and massless particles in Minkowski space [15] and on a 
general manifold [16]. These diffusions have an equilibrium limit determined by 
the friction. Linear diffusion equations are applied in the quark-gluon plasma 
and in heavy ion collisions [17] [18] [19] [20] (for non-linear diffusions see [21]). The 
relativistic invariance in these applications is not clear because only diffusion of 
some momentum variables is investigated (e.g. diffusion of the rapidity). The 
first diffusion equation applied to relativistic phenomena appeared in the paper 
of Kompaneets [5] (although the electron gas was treated in a non-rclativistic 
approximation and only the diffusion of the photon energy was considered). 
Non-relativistic non-linear diffusion equations appear in porous media, plasma 
physics, stellar dynamics and most importantly in hydrodynamics [22]. A non- 
linear equation in the phase space (Kramers equation) can be related to a non- 
linear equation in the configuration space in the high friction limit [23]. The 
statistical physics of non-linear diffusion equations has been discussed recently 
in [24] [25]. It seems that power-like non-linearities in the diffusion equation are 
intrinsically connected with the Tsallis statistics [26] [27] [28] [29] [30]. 

The organization of this paper is as follows. In sec. 2 we begin with rela- 
tivistic quantum mechanics in order to discuss the form of the nonlinear master 
equation. We treat Stiickelberg's proper time formalism in relativistic quantum 
mechanics [31] [32] [33] as a useful heuristic tool for a derivation of relativistic 
invariant equations. We consider the relativistic diffusion as a classical limit 
of the master equation. We propose a non-linear Lorentz covariant relativistic 
diffusion equation preserving the positivity and normalization of the probability 
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distribution. In sec. 3 we recall the model of the linear relativistic diffusion (we 
need this simple model to fix the notation). In sec. 4 we obtain a non-linear dif- 
fusive relativistic transport equation. We begin the studies of its consequences. 
First, we show that the transport equation has a stationary solution but no 
equilibrium. In sec. 5 relativistic covariant drifts are discussed. It is pointed out 
that an additional tensor is needed in order to define the drift. A unit four- 
vector which can be related to the frame velocity is introduced in order to 
construct the drift. The resulting diffusion equation has an equilibrium. A de- 
pendence of the physical equilibrium state on the relativistic frame of reference 
has been discussed in physics literature for a long time (see [34] and references 
quoted there). We consider quadratic and power- like non-linearities. In sec. 6 we 
show that the linear approximation coincides with the diffusion discussed in [15] 
and [35] (equilibrating to the Jiittner distribution, see also earlier papers cited 
in [7] [8]). In sec. 7 we define the relative entropy for the non-linear diffusions. 
We discuss drifts leading to the quantum equilibrium or to the Tsallis distribu- 
tion. We show that the relative entropy is decreasing monotonically when the 
diffusion tends to its equilibrium. 

2 Irreducible representations of the Poincare 
group and relativistic diffusion equations 

Relativistic invariance in quantum physics should be treated by means of the 
representation theory of the Poincare group (irreducible representations if ele- 
mentary particles are discussed). Let A € SL(2, C) and A be a homomorphism 
of SL(2, C) onto SO(3, 1). We consider one-particle states (described by a func- 
tion of the momentum) transforming under an irreducible unitary representation 
U of the Poincare group [36] [37]. When restricted to the Lorentz group the rep- 
resentation is expressed on eigenvectors |p, a) of the momentum P and the third 
component of the spin Sa(— 2j — 1 < <r < 2j + 1) 



where p^p^ = m 2 c 2 , 7Z is the Wigner rotation and the matrix V is defined by 
(2j + l)-dimensional representation of SU(2) . We obtain generators of the 
Lorentz group by differentiation (Greek indices denote space-time coordinates 
whereas the Latin indices refer to spatial coordinates) 



U{A)\p) = V(1Z(A))\A(A-')p), 



(1) 



(2) 



where 




(3) 



(4) 
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and Y,fj, v e su{2) . 

We outline here the scheme which leads to a relativistic diffusion equation 
discussed in more detail in [38] (the linear theory). For any classical observable 
O we can define its operator Weyl version O. Then, we consider matrix elements 
of O in the basis of the wave functions (1) 

0^(p,p') = (p,a\0\ P ',a'). 

We define a classical function on the phase space (the Wigner matrix) W as the 
Fourier transform of O (its limit h — > should coincide with O) 

W aa , (x, p) = / dkdk' J S(p - ±k - ±k') cxp(z(k - k')x)(9 ffff , (k, k'). (5) 

In the massless case instead of 2j + 1 states of spin j we have only two states 
of helicity A = j and A = — j. The formulae for generators of an irreducible 
representation of the Poincare group have been discussed in our earlier paper 
following [36] [39] [40] [41] [42]. 

We consider relativistic quantum mechanics formulated in an invariant way 
by means of the proper time r [31] [32] [33]. We look for no n- linear master equa- 
tions consistent with quantum mechanics. A non-linearity of the master equa- 
tion arises as a consequence of the quantum statistics [2] or in a description of 
particles interacting by means of an electromagnetic field [43] [44] [45] . In the 
first case, the probability of a transition depends on the occupation number of 
the state. In the second case an electromagnetic field depends on the distri- 
bution of charged particles which is described by the density matrix p. If we 
express the electromagnetic field by its source p then we obtain a non-linear 
master equation. In general, such an equation will have a complicated nonlocal 
form. We consider a local approximation independent of details of particular 
interactions. We assume that the operators in the equation can be expressed 
solely by generators of the Poincare group and some scalars. We demand that 
after a time evolution the density matrix remains a Hermitian operator and 
that its trace is preserved. We additionally assume that the generator of the 
evolution is a quadratic function of the generators of the Poincare group (this 
assumption leads to the diffusion). Then, the structure of the evolution equa- 
tion can be derived in a similar way as in the case of the Lindblad equation [46] . 
Our assumptions suggest the equation 

d T p = [P^p]] - y([M f " / , [GM^,p]} + [GM^, [M^p]]). (6) 

Here, j 2 is a diffusion constant which could be determined from the physical 
model of particle interactions, G(p) is an operator depending on p which is a 
scalar with respect to the transformations of the Lorentz group. Eq. (6) is 
covariant under transformations of the Lorentz group. If there is another tensor 
at our disposal then we may insert + in eq.(6). We shall discuss 
this possibility in sec. 5. 
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In the linear case (G = 1) the time evolution of an observable O can be 
defined by the time evolution of the state p (here the trace defines an expectation 
value of the observable) 

Tr(p T 0)=Tr{pO T ). (7) 

From the relation (7) between the evolution of states and the evolution of ob- 
servables it follows that an evolution equation for observables also preserves the 
positivity of an observable and its trace. Then, from eqs.(6) and (7) we obtain 
(see [38] and [3] [4]) a linear equation for W (defined in cq.(5)). We generalize 
this equation ( when G ^ 1) to a non-linear diffusion equation 

d T W = p»d*W - \j 2 (M^GM^W + WGM llv M ltv 
-GM^WM*" - M^WGM^) 

where G{W) is a Lorentz scalar but in general a non-linear function of W, the 
operators on the rhs of W act to the left (their position is relevant only for 
particles with spin, see [38]). 



3 Linear diffusion of particles without spin 

Let us write explicitly the linear version (G = 1) of eq.(8) for particles with the 
spin equal zero (see [38] for a higher spin) 

K- 2 (d T W - p»d*W) = |AgW, (9) 

where k 2 = "f 2 m~ 2 cr 2 and 

= m 2 c 2 (df +d 2 + fig) + p>p k d ] d k + ip k d k (10) 

k = 1,2,3 and dj = g^j. The derivatives over position will have an index x, 
derivatives without an index are over momenta. We can take the limit m — >• 
of the operator (10). Then, the limit m — > of eq.(9) is defined by 

A H =pjp k d j d k + 3p k d k . (ii) 

We can see that the operator (10) is the Laplace-Bcltrami operator on the 
hyperboloid p^p^ = m 2 c 2 

A™=g-id jg i k g^d k , (12) 

where 

g i k = m 2 c 2 5 jk +p>p k , (13) 

g = det(gj k ) and g^ 1 — m 6 c 6 p 2 ) , po = \J m 2 c 2 + p 2 . Note that in terms of the 
spatial components of the four-vector momentum the Lorentz transformation 
p — >• p' is non-linear. The invariance of the metric tensor takes the form 

9 jk (p') - Jjl9 lr (p)Jkr, (14) 
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where 

dp' 1 
~dp l 



J » = ±T- ( 15 ) 



Using eq.(14) we can check by direct calculations that transforming coordinates 
by means of the Lorentz transformation we do not change the Laplace-Beltrami 
operator (what is well-known; however, the explicit calculations checking the 
invariance of the non-linear diffusion will be useful later on) 

Ag(p') = AS(p). (16) 

In the classical description we define the evolution of the probability density 
Po& by an adjoint equation 

(* T , W) = / dxdp$ T (x,p)W(x,p) = J dxdp<S>(x,p)W T (x,p) = (*, W T ). 

(17) 

We define the adjoints of operators with respect to the (real) scalar product 
(17) 

($,AW) = (A*Q,W). (18) 

Note that L* k = —Ljk and LjJ • = idjPo- 

We obtain the transport equation if both sides of eq.(17) do not depend on 
t. This is the case if 

-P^W = yA£W (19) 

or equivalently 

p» d *$ = ^Ag**, (20) 

where A^ 1 * is the adjoint (18) of A^ (this is also the adjoint in L 2 (dp)). We 
write fi = po&. Then, eq.(20) can be rewritten in the form 

P^n = ^d j (^ k d k n). (21) 

4 Non-linear relativistic diffusion equations (no 
drift) 

We restrict ourselves to particles with spin zero from now on. Eq.(8) can be 
expressed by spatial components p of the four-vector p as 

8 n - Po 'pVxO - \n 2 d 3 (G{Q)g^Po 1 d k d^ (22) 

or in a more elegant way 

p^ = ^= g d j [^ k d k g) (23) 
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when 

G(fi) = g(ny. (24) 

Eq.(22) is in a divergence form 

A J A = 0, (25) 

where (x A ) = (x, p). As a consequence the particle number N is conserved 
(x = ct) 

d t N = d t J dxdpil = 0. (26) 

We could show that the preservation of the probability (26) is a classical equiv- 
alent of the quantum conservation laws 

d t Tr(p Q O t ) = d T Tr(0 T ) = 0. (27) 

We are looking for solutions of eq.(22) depending only on dimensionless energy 

e = /3cp , (28) 

where /3 _1 is a parameter of the dimension of the energy which we can set as 
equal to (fcsT) -1 (where T is the temperature and Ub the Boltzmann constant). 
Then, 

djfl(e) = f3cpjPQ 1 d e n(e). 
If 0, depends only on e then eq.(22) reads 

d t n t (e) = ( e - 2 9 £ (e 3 - em 2 c^ 2 )d t G(VL t ) + m 2 c 4 /3 2 e- 2 c> e S(ft t )) . (29) 

First, assume m > 0, then the stationary solution <9o^s = is determined by 
the equation (note that the current (25) J 7^ 0, hence fls is not an equilibrium) 

(e 3 - em 2 c i l3 2 )d f _g{^ s ) + m 2 c A p 2 g{n s ) = R, (30) 

where R is a constant. If Q is known then integrating eq.(30) we obtain an 
implicit equation for Sl s 

y/e z — m z c^fi z 

If the stationary state fis(e) is given then eq.(31) determines Q as a function of 
O. 

In the massless case, let the initial £1 depends on x and /3cp — ne. Then, 
eq.(22) reads 

9A(e) - cnV x O t = P^e- 2 d e e 3 d e g(n t ). (32) 



Let Ot(x, n, e) be the solution of eq.(29) (with m = 0, x and n treated as 
fixed parameters) then the solution of eq.(32) is fi t (x — net, n, e) . Solutions of 
eq.(30) have a discontinuity at m = 0. The limit m — > of eq.(31) exists but is 
trivial (Q = const corresponding to G = 0). There exist non-trivial stationary 
solutions of eq.(32) (which are not a limit of the solution (31)). In fact, in the 
massless case the condition d t flg = gives 

g(Q s ) = a- 2 e- 2 (33) 

where a is an arbitrary integration constant. 

Power-like nonlinearities in non-relativistic diffusion equations have been 
discussed in [26] [27]. A relativistic non-linear diffusion as a function of the 
rapidity has been studied in [21] and applied to the heavy ion collisions. In our 
model the Tsallis distribution [30] [47] 

n s = (l + (q-l)e)-^ (34) 
results as a stationary distribution if m = 0, a = q — 1 and 



Q = {q-iy 2 £l 2q - 2 (l-W- 1 )- 2 . (35) 

From the Lorentz invariance of the diffusion equation (22) it follows that if 
fis(e) is a stationary solution of eq.(22) then (Is(cj) is also a stationary solution 
of this equation where 

lu = Pcp^Wft (36) 

and w p is an arbitrary unit time-like four-vector. We can show a more general 
result: assume that ft depends only on u> and on X = x^w^ then eq.(23) can 
be expressed in the form 

d x n x (cj) = i^L^dj^ - ^^m 2 c 4 /3 2 )a w g(Ox) 

+w^w^m 2 c i l3 2 u- 2 d uj g{flx))- 

An A-independent solution of eq.(37) can be derived in the same way as in 
eqs.(30)-(31) with e — > u and m 2 c 4 — > m 2 c 4 w M w (1 .We recover eq.(29) when 
w = (1,0,0,0). 



5 Frame-dependent drifts of non-linear diffusion 
equations 

In sec. 2 we suggested that the non-linear relativistic diffusion equation is essen- 
tially unique up to an arbitrary (state dependent) diffusion strength G. In sec. 4 
we have shown that the relativistic invariant diffusion equation has a stationary 
solution Cls which is not an equilibrium (J ^ 0). In this section we construct 
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Lorentz covariant drifts such that the resulting diffusion equation has an equi- 
librium (in addition to the stationary states of sec. 4). We could begin (as in 
sec. 2) with the representation theory of the Poincare group. The first order 
operator transforms as a tensor. In order to define a drift we need another 
tensor to couple in an invariant way. As an example, an interaction with 
the electromagnetic field can be treated as a drift F^ v 'M M „W . However, if 
the equilibrium is to be space-time independent then such a drift is excluded. 
We cannot couple to M^ u for the same reason. The tensor p M p„ being sym- 
metric does not couple to M M!y . We need some other vectors (or tensors) to 
achieve such a coupling. In sec. 4 we have expressed the stationary state fls as 
a function of the energy. We have shown that a Lorentz transformation of fls 
depending on the velocity of the Lorentz frame is again a stationary solution 
of the diffusion equation. The velocity is determined in the unique way. We 
could assume from the beginning that is a variable in the diffusion equation 
(on physical grounds the equilibrium should depend on the Lorentz frame [34] ) . 
In such a case we obtain a wider class of covariant diffusion equations. We con- 
sider diffusion equations which in the rest frame and in a linear approximation 
reduce to the diffusion equations discussed in [ll][12][15](Boltzmann statistics). 
Assuming that there is an antisymmetric tensor A^ v available we suggest the 
following non-linear diffusion equation with a drift (M M „ = for a particle 
without spin, L*^ v is defined below eq.(18)) 

p^p^n = \nHL*^G(iLl v p^n + A„„(fi)po ^) • (38) 

In p coordinates eq.(38) reads 

y/gvr&pl = ±K 2 dj (Gg^ k ^g(d k n + A^Slfj . (39) 

We are looking for equilibrium solutions Q,e of eq.(39) depending on e. It is 
easy to see that such a solution is determined by the equation (the current J of 
cq.(25) is zero) 

g^p^dktlE + A°m E = (40) 

(the term A^ h is absent in eq.(40) because of the rotation invariance of SIe and 
the antisymmetry of A* k ). 

As an example of relativistic covariant eqs.(38)-(39) we could consider an 
interaction of a particle with an electromagnetic field (however, this does not 
serve our purpose of finding new equilibria) . In an electromagnetic field A^ v = 
Fpv = d^A v — d v A^ eq.(39) takes the form 

8 n - Po 'pVxfi = \^d-j (G W fc Po ^fi + F fc „p"po x 0) . (41) 

It can be shown that if fis(e) is a stationary solution of eq.(22) then the sta- 
tionary solution of eq.(41) in a scalar potential A is ^s(e + A ). 



9 



We could also approach the construction of a relativistic covariant drift di- 
rectly (without referring to the representation theory of the Poincare group). 
So, if is a scalar then it follows from eq.(14) that (for an arbitrary scalar 
H(Q)) the drift 

H(n)g jk d 3 (j)d k n 

is relativistic invariant. If the equilibrium is to be independent of x then we 
should construct the scalar from p and another vector (p 2 = m 2 c 2 is a 
constant). We put <p — w^p^ and identify u> M with the frame velocity. Such a 
choice is equivalent to setting 

A" v = p^p" - w v p ll )Q.- 1 H{Q) (42) 

in eqs.(38)-(39). Choosing a special frame w = (1, 0, 0, 0) gives A j — /3pjiJ(r2)f2 _1 
and Ajk = in eq.(39). Assuming that the equilibrium Qe (40) depends only 
on the energy (28) it is easy to see that there exists an equilibrium solution of 
eq.(40) if 

d e Sl E + H(Q E ) = 0. (43) 
Then, H(Q) = gives the Jiittner distribution [48] whereas 

H a (n) = 0(1 + an) (44) 

leads to the quantum distributions (er = +1 for bosons, a = — 1 for fcrmions 
and the Boltzmann statistics could be treated as a limit a — >• 0) . In general, 
the diffusion strength G could again be a function of O. 

In a rest frame, the diffusion equation (39) for scalar particles has the simple 
form 

d f2 - % 'pVxft = \^d 1 (G{^)(g^ k p^ 1 d k n + pcp>H{^y (45) 
In an arbitrary frame eq.(45) can be written in a covariant divergence form (25) 

VflP^^H = ^^(GHg^^^H^dkn + cM-pcp^dkexp^cp^w^y 

(46) 

If d^fls — in eq.(45) and J ^ then ils is called a stationary state. If 
J(Ob) = then Qe is the equilibrium. Eq.(46) has the equilibrium solution 
(this covariant form of the equilibrium probability distribution is discussed in 
particle physics in [49], see also [34]) 

n E = (zexpiPcufpn) - <t) 1 (47) 

for H defined in eq.(44)( here z = exp(— /x) and [i is the chemical potential; the 
Jiittner distribution corresponds to a = 0). 
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The equation with the drift H q (Sl) = SI 9 

d Q n - Pq x pV x n = \K 2 djG{Q) (grtp^dktt + /?cpj-fi«) (48) 

gives the Tsallis equilibrium distribution (34). 

We would like to find a quantum analog of the Tsallis distribution. There are 
already some candidates in the literature [50] . We suggest another one resulting 
from the drift 

H aq = + aSl)y ({1 + afl) q - aSl^ 1 (49) 
Eq.(43) gives an implicit equation for the equilibrium 

n x - q -{l + aSl E ) 1 - q = («-l)(e-/i). (50) 

Eqs.(49)-(50) have the correct limit q — > 1 (quantum distribution) and a — > 
(classical distribution). We show in sec. 7 that the drift H aq allows to define a 
monotonic relative entropy. 

If we assume that ft depends only on e then eq.(45) takes the simple form 



flA(e) = ^e- 2 (d £ ((e 3 - m 2 c 4 /3 2 e)c> e S(Q t )) + 
m 2 c 4 /3 2 9 e £(Q t ) + m 2 c 4 p 2 GH + 5 £ ((e 3 - m 2 c 4 ,3 2 e)Gtf 



(51) 



The stationary distribution is determined by the equation 9 t f2 = 0. After an 
integration of the rhs of eq.(51) we obtain the first order equation for Slg 

where R is an arbitrary integration constant. R = gives the equilibrium (43). 
If H = SI (corresponding to the Juttner distribution ) and G = 1 then eq.(52) 
is linear. It can be checked that the only normalizable (J Q$ < oo) solution is 
at R = 0. However, there may be normalizable solutions of eq.(52) with R ^ 
if G 7^ const. The obvious example is G ~ 

In the limit m — > eq.(45) takes the form resembling the Kompaneets equa- 
tion [5] 

dtQ - cnV x fi = ^-e- 2 -^- (e 3 G(r>) (^ft + if (ft)) ) . (53) 

where n = p|p| _1 . 

The stationary distribution is determined by 

-n s + H(Q s ) =Re- 3 G(Sl s )- 1 . (54) 

The Kompaneets equation has been derived in quantum electrodynamics of 
clcctron-photon scattering [5] [6] [51] with a non-relativistic approximation for the 
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electron velocity distribution. In such factor e 4 enters eq. (53) instead of 

e 3 . In refs.[52][53] e 4 is replaced by a general function a(e) in order to take rel- 
ativistic corrections into account. Our discussion based on eqs.(14)(39)and (46) 
shows that the complete differential equation in the phase space is relativistic 
covariant only when a(e) — e 3 . The Kompaneets equation has important appli- 
cations in astrophysics and cosmology (Sunyaeev-Zeldovitch effect). It is also 
studied as a mathematical model for the time dependence of the Bose-Einstein 
condensation [51] [54] [55]. 

6 Low density approximation (linear diffusion) 

We recall the linear diffusion of ref.[15] and its thermodynamics [35] in order to 
compare it to the non-linear one. If we assume that f2 is small then the general 
drift H can be expanded in a Taylor series 

H(si) = n + a 2 n 2 + ... 

(H(0) = if the solution of eq.(43) is to have the required property Q,e — > 
for e — > oo ). We assume in this section that G(SY) ~ 1 and f2 >> il 2 . Then, 
eq.(45) becomes linear 

d Q n - Pq x pV x O = i^ 2 ^ (g> k Po 1 d k Sl + PcpjSl) . (55) 

It has the Juttner equilibrium solution [48]. 

We wish to describe the time evolution in terms of some extensive (ther- 
modynamic) functions of the probability density. It seems that the relative 
entropy measuring the distance between two probability distributions plays the 
crucial role in such a description (the relative entropy measures the speed of 
convergence f2 t — > Qe, its relevance has been discovered in [56] [57]; the relative 
entropy in nonextensive statistical mechanics is discussed in [30] ) We define the 
Boltzmann-Kullback-Leibler relative entropy Sk of two unnormalizcd probabil- 
ity distributions Q and Qe (N and Ne are the normalization constants) 

s K {a- t n B ) = n- 1 J dxdpmn^- 1 ^^)- 1 ^). (56) 

It is known that [58] 

s K (n ; n E )>o. (57) 

An easy calculation using the diffusion equation (55) gives 

d s K {n: q e ) = -\n- x j dpdxpo 1 n- 1 gi k (d j n + pc PjP ^n){d k n + pep^n). 

(58) 

It follows from eqs.(57)-(58) that Sk{^',^e) is a non-negative function mono- 
tonically decreasing to zero when t — > oo. 
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We define the Boltzmann entropy 

S(Q) = -k B j dxdpfilnfi (59) 

the energy 

£ = c j dpdxpofl (60) 

and the free energy 

T = r 1 N(s K (n ; n E )-\n(NN E : 1 )). (61) 
Then, from eqs.(58)-(61) we obtain the basic thermodynamic relation 

TS = £-T. (62) 

7 Thermodynamics of the non-linear diffusion 

In general, a meaningful (non-Hamiltonian) evolution equation should lead to 
an increase of entropy. Conversely, if the entropy is defined then the evolution 
equation could be determined from the requirements that the entropy is in- 
creasing and bounded. If we define the entropy S in terms of an entropy density 
s 

S = - J cixdps(O) (63) 
then for the diffusion (22) (without drift) 

d ° S ^l J dxd PPa ls "(^)G(^)d^9 jk dk^- (64) 
It follows that S is increasing if 

s" > 0. (65) 

The diffusion without a drift has no equilibrium but only a stationary state. 
There may be no limit for the increase of the entropy. For diffusions with the 
drift (44) quadratic in f2 (quantum equilibrium distributions) as well as for 
H q = fi 9 and H aq (49) (Tsallis distributions) we can define the relative entropy 
Sk{&, &e) which is monotonically decreasing to the equilibrium state when 
t — > oo (Sk of eq.(58) would not be a monotonic function). For the drift (44) 
define 

S^(Q,n E ) = J dxdp^ln^r^ 1 ) - a(l + cxO) In ((1 + aO)(l + an E )-^ 

(66) 
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for H = Sl q (such a relative entropy is suggested also in the Appendix of 
for a general theory see [30]) 

s 9 K (si,si E ) = (2 - q)-\l - q)- 1 J dxdpn 2 -i 
-(i - q)- 1 J dxd P nn 1 E q + (2 - q)- 1 J dxdpQ 2 E q . 

We define Tsallis relative quantum entropy for the drift H aq (49) 



(67) 



s a K q {fi,n E ) = (2-q)- 1 (i- q )- 1 J dxd P n 2 -i 

-(1 - q)- 1 J dxdpfin 1 E q + (2 - q)- 1 J dxdpSl 2 E q 

-a{2 - q)' 1 ^ - q)- 1 f dxdp(l + aSl) 2 - ( ' + a(l - q)' 1 j dxdp(l + aSl){l + aSl E ) 1 - ( i 
-a(2 - q)- 1 J dxdp(l + an E ) 2 - q . 

(68) 

The three entropies (66)-(68) have the properties Sk(SI E , Sl E ) = , S' K (Sl E , Sl E ) = 
and S K (Sl,Sl E ) > 0. It follows that Sk(SI, ^e) > 0(in contradistinction to 
Sk in eq.(58) we do not need to normalize SI and Sl E in order to prove the 
positivity of Sk in (66)-(68)) . 
In an equilibrium state 

r = G(si E ) g rk Po 1 (d k si E + f3cp kP ^H(si E )) = o. (69) 

Then, using eq.(39) for the calculation of the time derivative of Sk we obtain 

d s K (si,si E ) = -U dxdvpu 1 H{£i)- 1 G(n) g 3 k {d ] n + A mP »n)(d k ti + A k ^n)) 

(70) 

(the Lorentz invariance of eq.(70) follows from eq.(14)). In the rest frame 

d s K (si, si E ) = -U d X dp P o 1 h{si)- 1 g{si) 



^' fc (a,n + /3cp,Po r H(fi))(9 fe O + pcp k p^H{Sl))) 



where H is either H a ,H q or if^. 



In the massless case eq.(71) simplifies to 

d s K {si,si E ) = -\ J dpp^HW^GffltfdjSi + pc Pa H{n)) 2 . (72) 

It follows that Sk{SI, Sl E ) > and the time evolution will decrease Sk until 
it achieves its minimum at the equilibrium (then Sk(SI E ,SI e ) = 0). If the 
state depends solely on the energy then ( this equality has been derived in the 
particular case of the Kompaneets equation in ref.[53]) 

d S K (Q, Sl E ) = -4tt J dee 3 GH- 1 (d £ n + H) 2 . (73) 

In the stationary state Sis (54) (when J ^ 0) there is a correction to the formulas 
(70)-(73). So, in the massless case under the assumption that both SI and Sis 
depend only on e we obtain 

doS K (Sl,Sls) =-4tt / dee 3 G(Sl)H{Sl)- 1 (d c Sl + H{Sl)) 2 

+4ttRJ deG^GiSls)- 1 HiSlsy^d.Sl + H(Sl)). ( ' 
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The i?-correction has an indefinite sign. Eq.(74) shows that Sk(^1,Qs) is not 
monotonic. fit does not tend to fig. It is convergent to SIe- 

We define the entropy (as suggested by eqs.(66)-(68)) for the Bose-Einstein 
and Fcrmi-Dirac distributions 

Sv^-keJ dxdp(n\nn-a(l + an)\n(l + <7n)^ (75) 

and for the quantum Tsallis distribution 

S aq = -k B (2 - q)-\l - q)- 1 J dxdptf-i 

+k B a{2 - q)- 1 ^ - q)- 1 J dxdp{l + af!) 2 ^ { '°> 

(when a = we have S aq — > S q ). In the models discussed in this section the 
principle of the maximum of entropy is satisfied. S is achieving the maximum 
(with fixed JO = 1 and J efl = £) at the equilibrium state He- We define the 
free energy as T = £ — TS ( where £ is defined in eq.(60)). From eqs.(66)-(68) 
and (75)- (76) we obtain that T = Sk + Te where Te is the time independent 
equilibrium free energy. Then, eq.(62) is satisfied as an identity. 



8 Discussion and summary 

The non-linear diffusion equations result from a macroscopic averaging of com- 
plex microscopic phenomena. We have shown that the basic requirements im- 
posed on such equations determine their form. The linear part of the equation 
is already determined by the Lorentz invariance and the requirement of the pos- 
itivity of the probability density. The non-linearity of the equation comes from 
quantum statistics (dependence of the transition rate to a state on the occu- 
pation of the state) and from an elimination of some fields which are averaged 
in the reduced dynamics. We have shown that the form of the non-linearity 
is closely related with the formula for the entropy. It is possible to obtain a 
drift for a linear diffusion such that the diffusion is equilibrating to the Bose- 
Einstein or Fcrmi-Dirac equilibrium distribution [15]. However, it seems that 
the non-linear drift is indispensable if we wish to define the relative entropy as 
a monotonically decreasing function which satisfies the thermodynamic relation 
(62). We can apply the non- linear diffusion equation to a statistical description 
of the motion of a stream of relativistic particles in a medium of some other 
relativistic particles. Such a problem is studied in high-energy physics [3] and 
astrophysics [60] [61]. However, the whole space-time dependence of the phase- 
space distribution is not discussed in these papers. It can be a complicated 
problem to derive an equation of space-time evolution of the stream of particles 
on the basis of a theory of fundamental interactions (see [6] [62] [63] [64] ) . It may 
require non-perturbative methods. In general, the resulting equation will be 
non-local. If the Markovian approximation is applicable to the full space-time 
dynamics then the relativistic invariance leads to the equations discussed in this 
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paper. The result of the multi-particle ultra-relativistic scattering may be in- 
dependent of the details of the interaction but depend more on kinematics and 
the relativistic processes of dissipation. In such a case we could check in ex- 
periments assumptions underlying the dissipative dynamics of cq.(6). Diffusion 
equations already served as a theoretical basis for RHIC data analysis [20] [21]. 
The diffusion model could be tested in heavy ion collision experiments. The 
complete space-time diffusion equation is also interesting as a model for space- 
time evolution of the Bose-Einstein condensation (discussed in atomic physics 
[54] and in astrophysics as a model for the star formation from the dark matter 
[65]). 
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